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15 Abstract 

16 Coronaviruses from both the Alpha and Betacoronavirus genera, interfere with the type I interferon 

17 (IFN) response in various ways, ensuring limited activation of the IFN response in most cell types. 

18 Of Gammacoronaviruses that mainly infect birds, little is known about activation of the host 

19 immune response. We show that the prototypical Gammacoronavirus, infectious bronchitis virus 

20 (IBV), induces a delayed activation of the IFN response in primary renal cells, tracheal epithelial 

21 cells and in a chicken cell line. Ifnfi expression in fact, is delayed with respect to the peak of viral 

22 replication and accompanying accumulation of dsRNA. In addition, we demonstrate that MDA5 is 

23 the primary sensor for Gammacoronavirus infections in chicken cells. Furthermore, we provide 

24 evidence that accessory proteins 3a and 3b of IBV modulate the IFN response at the transcriptional 

25 and translational level. Finally, we show that, despite the lack of activation of the IFN response 

26 during the early phase of IBV infection, signalling of non-self dsRNA through both MDA5 and TLR3 

27 remains intact in IBV-infected cells. Taken together, this study provides the first comprehensive 

28 analysis of host-virus interactions of a Gammacoronavirus with avian innate immune responses. 

29 

30 Importance: Our results demonstrate that IBV has evolved multiple strategies to avoid activation 

31 of the type I interferon response. Taken together, the present study closes a gap in the 

32 understanding of host-IBV interaction, and paves the way for further characterization of the 

33 mechanisms underlying immune evasion strategies as well as pathogenesis of 

34 Gammacoronaviruses. 
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Introduction 

Coronaviruses constitute a large family of positive-stranded RNA viruses and cause a range of 
human and veterinary diseases. Infectious bronchitis virus (IBV) is the prototype avian coronavirus 
from the Gammacoronavirus genus and the causative agent of a highly contagious respiratory 
disease of major economic importance to the poultry industry (1). IBV enters the avian host 
through the respiratory tract, where it causes destruction of the epithelium leading to respiratory 
distress and initiation of secondary bacterial infections. Depending on the strain, IBV can also 
spread to other epithelial surfaces such as the gastrointestinal tract, the kidneys and the oviduct, 
the latter causing problems in egg production and quality (1-6). Contrary to coronaviruses from the 
Alpha and Beta genera, including human coronavirus HCoV-229E, Severe Acute Respiratory 
Syndrome (SARS-CoV), Middle East Respiratory Syndrome (MERS-CoV) and mouse hepatitis virus 
(MHV), very little is known about how Gammacoronaviruses including IBV evade or interfere with 
innate immune responses of their host. 

Innate immune responses consist of a network of antimicrobial mechanisms, of which the type I 
interferon (IFN) response is an essential defence mechanism against viruses. Typically, the type I 
IFN response, from hereafter referred to as IFN response, is initiated upon activation of host 
pattern recognition receptors (PRRs), present in all animal cells. Two families of PRRs have been 
shown to be involved in the recognition of RNA viruses namely the membrane-bound Toll-like 
receptors (TLRs) and the cytosolic RIG-I-like receptors (RLRs) (7). The primary ligands for 
activation of these PRRs are double-stranded RNA (dsRNA) and 5' triphosphate-containing RNA, 
normally absent from uninfected host-cells. Activation of RLRs leads to the transcription of genes 
encoding type I interferons (IFNa and IFN|3). These interferons are secreted from the infected cell 
providing a signal for the infected as well as the neighbouring cells that induce the transcription of 
anti-viral effector genes collectively called interferon stimulated genes (ISGs). 
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The ability of a virus to replicate and produce infectious progeny depends for a large part on its 
ability to avoid induction or counteract the IFN response of its host. Indeed, a common feature of 
Alpha- and Betacoronaviruses, including FICoV-229E, SARS-CoV, and MFIV, is their limited 
activation of the IFN response (8-13). This limited activation can be partially explained by 
intracellular membrane rearrangements that might shield dsRNA and other viral components from 
recognition by host PRRs (14, 15). In addition, coronavirus nspl6 displays 2'-0-methylase activity, 
which results in 2'-0-methylation of a ribose moiety on the 5' cap of coronavirus mRNAs, making 
them indistinguishable from host mRNAs (16). Furthermore, many other coronavirus proteins, such 
as nspl, nsp3, the nucleocapsid and many of the accessory proteins have been shown to interfere 
with the IFN response in various ways (reviewed in (17, 18)). 

Interaction between G ammacoronaviruses and innate immune responses of their avian hosts is 
poorly understood. Early studies on Gammacoronaviruses in chicken suggest that IBV-induced IFN 
production is variable and dependent on both virus strain and cell type. (19-22). Further, two 
transcriptional studies on tissues collected after in vivo and in ovo IBV infections, found only limited 
upregulation of ISGs at 1 - 3 days post-infection (23-25). Functional studies using IBV Beaudette 
showed that it induced cell-cycle arrest and apoptosis (26, 27), that IBV interacts with eIF3f (28) 
and that IBV inhibits protein kinase R activation, thereby maintaining protein synthesis (29). 
Although these studies did provide a number of details on the interactions between IBV and the 
host cell, most experiments were carried out in Vero cells. This non-avian cell line is one of the 
very few cell lines in which the IBV-Beaudette strain has been adapted to grow, facilitating in vitro 
experiments. Vero cells, however, lack the Ifn/3 gene, preventing them from mounting a type I IFN 
response (30, 31), reducing the value of Vero cells for research on innate immune responses to 
IBV. In addition, the Beaudette strain is non-pathogenic in vivo with limited replication in host 
tissues (32), reducing the value of these in vitro studies for translation to in vivo situations. For 
these reasons, we used pathogenic isolates of IBV to infect primary chicken cells, and a chicken cell 
line, as these isolates are known to infect, spread and cause clinical disease in vivo. 

In the current study, we show that IBV infection leads to a significant induction of Ifn/3 transcription 
through an MDA5-dependent activation of the IFN response, albeit delayed with respect to both 
virus replication and accumulation of dsRNA. This delayed induction of Ifn/3 was further confirmed 
through RNA FISH analysis showing that accumulation of Ifn/3 mRNA is restricted to IBV-infected 
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and not neighbouring uninfected cells. Although the time lag between accumulation of dsRNA and 
induction of Ifn/3 transcription might suggest that IBV interferes with recognition of dsRNA, we 
observed that sensing of exogenous (non-self) dsRNA remained functional in IBV-infected cells. 
Using mutant IBV viruses we demonstrate that both accessory proteins 3a and 3b are involved in 
limiting Ifn/3 expression, as both 3a and 3b null viruses induced increased Ifn/3 expression. 
Nevertheless, 3a and 3b seem to have a differential effect on IFN protein production, infection with 
3a null virus induced lower IFN levels whereas a 3b null virus increased IFN production compared 
to the parental virus. Altogether, our data suggest that IBV delays but does not prevent detection 
by MDA5, and that accessory proteins 3a and 3b modulate the IFN response in avian cells. This is 
the first study addressing immune evasion and interference strategies of IBV in chicken and not in 
mammalian cells, providing information essential to further understanding of the pathogenesis of 
Gammacorona viruses. 

MATERIALS AND METHODS 

Cells 

Chicken embryonic kidneys were aseptically removed from 17- to 19-day-old chicken embryo's 
(Charles River, SPAFAS). A cell suspension was obtained by trypsinisation for 30 min at 37°C and 
filtered through a 100 pm mesh. The resulting chicken embryo kidney (CEK) cells were seeded at 4 
x 10 5 cells/cm 2 in 199 medium (Invitrogen) supplemented with 0.5% fetal bovine serum (FBS, 
FASC) and 1% PenStrep (Gibco®, Invitrogen). Chicken trachea cells were isolated from 8- to 10- 
week-old chickens (white leghorn). Tracheas were collected in ice-cold PBS, washed and stripped 
from adipose tissue. Trachea were filled with a solution of 3.5 U/ml protease type XIV (Sigma), 4 
U/ml DNase I (Qiagen) and 1% PenStrep in EMEM, sealed with clamps, and incubated overnight at 
4°C. The next day, cells lining the luminal side of the trachea were flushed out with cold EMEM, 
filtered through a cell strainer and seeded at 4 x 10 5 cells/cm 2 in DMEM supplemented with 10% 
FBS and 1% PenStrep. The RIG-I wt , RIG-I K0 , MDA5 wt and MDA5 K0 MEFs were provided by Prof. S. 
Akira (33). The MAVS wt and MAVS K0 MEFs were provided by Z.J. Chen (34). DF-1, CEC-32 and 
MEFs cells were cultured in DMEM (Gibco®, Invitrogen) supplemented with 10% FBS and 1% 
PenStrep. All cells were incubated in a humidified incubator at 37°C and 5% C0 2 . 
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DF-1 Ifnp-luc reporter cell line 

DF-1 cells were transfected using Fugene (Promega) according to the manufacturer's instructions 
with a construct expressing firefly luciferase under the control of the -llObp proximal region of the 
human 7F/V/3 promotor (35). Stably expressing cells were selected over a period of 3 weeks using 
geneticin (500 pg/ml). DF-1 Ifnp-luc stable cells were cultured in DMEM supplemented with 10% 
FBS and 1% PenStrep and were not further subcloned. 

Viruses 

IBV-M41, IBV-QX and IBV-Italy-02, Rift Valley Fever Virus clone 13 (RVFV CI13) and Infectious 
Pancreatic Necrosis Virus (IPNV) were obtained from Merck Animal Health, Boxmeer, The 
Netherlands. Sindbis-GFP was a kind gift from J. Fros, (Laboratory of Virology, Wageningen 
University). IBV Beaudette, strain Beau-R, as well as the generation of the ScAUG3a, ScAUG3b, 
ScAUG3ab and ScAUG5ab Beau-R null viruses has been published previously (36-38). In these 
mutant IBV viruses, the start codons of the indicated accessory genes were mutated to stop 
codons. All IBV strains were amplified and titrated on CEK cells. Sindbis-GFP was amplified on baby 
hamster kidney (BHK) cells and titrated on CEK cells. RVFV Cl 13 was amplified and titrated on Vero 
cells, an African green monkey cell line. IPNV was amplified and titrated on the CHSE-214, 
Chinook-salmon cell line. IPNV was inactivated by 20 min UV exposure on a 48W BXT-26-M 
(Uvitec). 

Poly I:C stimulation and RNase treatment 

Polyinosinic-poly(C) [p(I:C)] sodium was purchased from Sigma, dissolved in nuclease-free water 
and stored at -80°C. p(I:C) was either directly added to the medium or transfected using 
Lipofectamine 2000 (Invitrogen) as per manufacturer's instructions. DF-1 cells (3 x 10 5 /well) were 
cultured in 24 well plates and transfected with 500 ng p(I:C). RNase treatment of CEK cell culture 
supernatant was performed by addition of 10 pg/ml RNase A (Invitrogen) before IBV infection or 
before stimulation with 2 pg/ml p(I:C). 

RNA isolation and cDNA synthesis 

Approximately 8 x 10 5 CEK cells or 3 x 10 5 DF-1 cells were lysed in RLT buffer (Qiagen) at various 
time points after treatment or infection. RLT cell lysis buffer was spiked with 1 ng/sample of 
luciferase mRNA (Promega) immediately prior to RNA isolation. Luciferase expression will later be 
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used as external reference gene for normalization during the gene expression analysis. Total RNA 
was isolated using the RNeasy Mini Kit (Qiagen) according to the manufacturer's instructions, 
including an on-column DNase treatment with RNase-free DNase (Qiagen). Before cDNA synthesis 
of 0.5-1 pg total RNA, a second DNase treatment was performed using DNase I, amplification 
grade (Invitrogen). Synthesis of cDNA was performed using Superscript III (Invitrogen) using 
random primers. cDNA samples were further diluted 1:50 in nuclease-free water before real-time 
quantitative PCR analysis. 

Gene expression analysis 

Real-time quantitative PCR was performed on a Rotor-Gene 6000 (Corbett Research), using 
Brilliant SYBR Green quantitative PCR (Stratagene) and primers (39-42) as listed in Table 1. Cycle 
thresholds and amplification efficiencies were calculated by the Rotor-Gene software (version 1.7). 
The relative expression ratio of the target gene was calculated using the average reaction efficiency 
for each primer set and the cycle threshold (C t ) deviation of sample vs. control at time point Oh, as 
described in (43). For calculation of the fold change of IBV total RNA, C t deviation was calculated 
versus C t 30, as no IBV was present in the non-infected cells that were used as control in all the 
experiments. Because expression of various housekeeping genes was unstable during virus 
infections at time points later than 24 hours (data not shown), gene-expression ratios were 
normalised using an external reference gene (luciferase). 

Immunohistochemistry 

CEK cells were seeded on fibronectin-coated glass Biocoat coverslips (BD Biosciences) at a density 
of 1 x 10 5 cells/cm 2 . After incubation at 37°C for 48 hours, cells were infected with IBV strain M41 
at an MOI of 1, and fixed at different time points with 3.7% paraformaldehyde and permeabilized 
using 0.1% Triton X-100. Infected cells were probed with anti-dsRNA antibody (English & Scientific 
Consulting) and polyclonal chicken serum raised against IBV M41 was obtained from Merck AH. 
Detection was done performed using Aiexa 488 goat anti-mouse antibody (Invitrogen) and FITC 
labelled goat anti-chicken antibody (Kirkegaard and Perry laboratories). Nuclei were stained with 
4',6-diamidino-2-phenylindole (DAPI). Cells were imaged using a Zeiss Primo Vert microscope and 
Axiovision software. Image overlays were made in ImageJ. 
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RNA fluorescence in situ hybridization 

RNA fluorescence in situ hybridization (FISH) was performed according to previously described 
protocols (44-46). A set of forty RNA FISH probes (20 bp), each labelled with one CAL Fluor® Red 
610 fluorophore and targeting chicken Ifni 3 (ENSGALT39477), was designed using the Stellaris® 
probe designer (Biosearch Technologies; https://www.biosearchtech.com/stellarisdesigner/). The 
coding sequence of chicken Ifn(3 is 601 bp, therefore to accommodate the optimum number of 
fluorescent probes (48, as explained in reference (44)), the 3'UTR was included in the probe design 
tool. CEK cells were grown on fibronectin-coated coverslips (BD Biosciences) at a density of 2 x 10 5 
cells/cm 2 . After incubation at 37°C for 48 hours, cells were infected with IBV M41, and at the 
indicated time points fixed in 70% ethanol at 4 °C. Hybridisation of the probes was performed 
using the manufacturer's protocol for adherent cells. Imaging was performed using a Roper (Evry, 
France) Spinning Disc Confocal System on a Nikon Eclipse Ti microscope using a 100 x Plan apo oil 
immersion objective (NA 1.4) and a 491 nm laser line. Z-stacks were collected with 0.25 pm Z- 
intervals. For each channel, maximum Z-stack projections were made and processed with ImageJ. 

Chicken type I IFN bioassay. 

Bioactive chicken type I interferon (chIFN) was measured using a bioassay based on the CEC-32 
quail reporter cell line expressing luciferase under the control of the chicken Mx promotor (47) 
(kindly provided by Prof. Peter Staeheli). Briefly, CEC-32 were incubated with serial dilutions of 
chIFN-containing samples for 6 hours, after which luciferase activity was quantified and IFN 
concentrations calculated using a chIFN standard. To avoid influence of IBV on the assay, samples 
were heat inactivated at 56°C for 30 min, which did not influence type I chIFN bioactivity. 

Gene silencing 

siRNAs targeting chicken Tlr3 and Mda5 were designed by and purchased from Microsynth, 
Switzerland (sequences in Table 1). Transfections were performed using siLentFect (Biorad) at a 
final siRNA concentration of 20 nM. For one well, 160 ng siRNA was combined with 1 pL siLentFect 
in 100 pL OptiMEM (Gibco) and incubated for 20 min. The siRNA complexes were added to 2 x 10 5 
DF-1 cells grown in 500 pL medium per well in a 24 well plate. siRNA complexes were left on the 
cells for 48 hours before further experiments were performed. 


226 Statistics 

227 All statistical analyses were performed in GraphPad Prism 5.0. RT-qPCR fold changes were first log 

228 transformed and then used for statistical analysis. For all tests, equality of variance was assessed 

229 using Bartlett's test. Significant differences (P<0.01) were determined by a one way or two-way 

230 ANOVA (indicated in the figure legend) followed by a Bonferroni post-hoc test. 
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Results 


IBV delays the onset of an IFN response during infection of primary chicken cells 

To investigate the kinetics of viral replication and IFN induction upon infection with the avian 
Gammacoronavirus IBV, we infected primary CEK cells (24) with the IBV M41 strain. To monitor 
the kinetics of the IFN response in relation to IBV replication, we quantified transcription of Ifn/3, a 
set of genes involved in innate immunity, extracellular IFN protein production, virus titres and IBV 
RNA in M41-infected CEK cells. In line with previous observations (48), progeny virus was produced 
after 6 hpi and virus titres reached a maximum around 24 hpi (Fig. 1A). Total intracellular IBV RNA 
levels reflected the kinetics of infectious IBV virus in the supernatant (Fig. 1A), reaching maximum 
levels around 24 hpi. Ifn/3 expression was delayed with respect to the peak of viral replication and 
remained low until 18 hpi, after which it was strongly upregulated, peaking around 36 hpi (Fig. 
IB). IFN protein activity levels were quantified using a chicken IFN-specific Mx-luc cell-based 
bioassay showing accumulation of IFN from 36 hpi onwards (Fig. IB). Concomitant with Ifn/3, a 
subset of genes involved in innate immunity, including Mx, Oas and 118, were upregulated whereas 
others, such as Tlr7, Adar, Isgl2, MHC-I and Ifnar2 appeared not, or only marginally affected by 
IBV infection (Fig. 1C). Pattern recognition receptors Mda5 and Tlr3 and the transcription factor 
Irf3 were also upregulated (Fig. 1C), which is of interest given the role of these PRRs in virus 
recognition. 

The delayed IFN response is independent from the cell type or virus strain 

Ifn/3 transcription during infection with coronaviruses such as MFiV and SARS-CoV is generally low 
(9, 10, 12, 13, 49), and was shown to be dependent on cell type and virus strain (50). The delayed 
induction of Ifn/3 transcription observed in IBV M41-infected CEK cells prompted us to investigate 
whether induction of Ifn/3 would be dependent on the cell type or IBV strain. Epithelial cells isolated 
from trachea of 10-week-old SPF chickens and DF-1 chicken fibroblast cells were infected with IBV 
M41 or IBV Beaudette (Beau-R, (38)). At several time points after infection, Ifn/3 levels were 
monitored by RT-qPCR (Fig. 2A and 2B). In both cell types Ifn/3 transcription followed the same 
kinetics observed in CEK cells (Fig. IB), indicating that induction of Ifn/3 by IBV is independent of 
cell type. To study whether induction of Ifn/3 transcription differs between different strains of IBV, 
we also infected CEK cells with the QX and It02 strains of IBV (Fig. 2C). Although we observed 
some differences in absolute levels of Ifn/3 upregulation induced by QX, It02 and M41, kinetics of 
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Ifn/3 transcription were similar, suggesting that delayed induction of Ifn/3 transcription could be 
considered a general feature of IBV infection in chicken cells. 

To assess whether CEK and DF-1 cells do have the intrinsic ability to express Ifn/3 earlier than 18h, 
we stimulated these cells with extracellular polyI:C (pi:C), transfected pI:C (t[pl:C]) or with the 
dsRNA virus Infectious Pancreatic Necrosis Virus (IPNV). We found that stimulation of CEK cells 
with pI:C could induce Ifn/3 transcription as early as one hour after stimulation (Fig. 3A). In DF-1 
cells, stimulation with IPNV, and t[pI:C], but not pI:C, induced Ifn/3 already at 4h (Fig. 3B). The 
observation that DF-1 cells do not respond to stimulation with extracellular dsRNA, is in accordance 
with previous findings and is most likely due to the lack of surface expression of TLR3 (51). In 
addition, a 12h infection of CEK cells with Sindbis, IPNV or Rift Valley Fever Virus clone 13 (RVFV 
Cl 13) induced a clear transcription of Ifn/3 (Fig. 3C). These results suggest that delayed expression 
of Ifn/3 is a specific feature of IBV infection and not an intrinsic characteristic of chicken cells. 

The intracellular pattern recognition receptor MDA5 is the primary sensor of IBV 

In general, dsRNA has been shown to be the canonical inducer of Ifn/3 during infection with Alpha- 
and Betacoronaviruses (16, 50). To determine which pattern recognition receptor (PRR) would be 
involved in sensing (ds)RNA of the Gammacoronavirus IBV, leading to subsequent Ifn/3 
transcription, we first examined the possibility that IBV-(ds)RNA could be sensed extracellularly by, 
for example, cell-surface receptors. To investigate this, CEK cells were infected with IBV M41 in the 
presence of RNase A and Ifn/3 expression was analysed. As a positive control, CEK cells were 
stimulated with pI:C in the presence or absence of RNase A. The IFN response to pI:C was greatly 
inhibited by addition of RNase A, which had no effect on Ifn/3 levels induced by infection with IBV 
M41 (Fig. 4A). These data suggest that Ifn/3 upregulation during the late stage (>18 hpi) of IBV 
infection could be the result of sensing of IBV-(ds)RNA by an intracellular rather than an 
extracellular pattern recognition receptor. This is consistent with our observation that IBV infection 
can be detected by DF-1 cells, which show only a marginal upregulation of Ifn/3 transcription in 
response to extracellular dsRNA (see Fig. 3B). In general, dsRNA can be recognised by membrane- 
bound TLR3 and cytosolic RLRs such as MDA5 and RIG-I. Genome mining strongly indicates that 
chickens do not express a RIG-I homologue (52), leaving TLR3 and MDA5 as the two PRRs 
potentially involved in dsRNA sensing. Silencing of Mda5, but not Tlr3, in DF-1 cells resulted in a 
70% decrease in Ifn/3 transcription (Fig. 4B). Similar results were obtained with an Ifn/3- luc DF-1 
reporter cell line in which silencing of Mda5, but not Tlr3, resulted in a 70% decrease in luciferase 
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activity by the reporter cells (Fig. 4C). Because no antibody against chicken MDA5 is currently 
available for protein detection, successful knockdown was evaluated using RT-qPCR demonstrating 
a silencing efficiency for both Tlr3 and Mda5 of approximately 60% (data not shown). Replication of 
IBV at the investigated time point was not affected by knockdown of neither Tlr3 nor Mda5, as 
measured by both virus titre and intracellular IBV total RNA (Fig. 4D). These results indicate that 
MDA5 is the primary PRR responsible for sensing Gammacoronavirus IBV-(ds)RNA in chicken cells. 

Early accumulation of dsRNA in IBV-infected cells does not result in early induction of 

Ifnp 

Having assessed that chicken cells can indeed promptly respond to stimulation with dsRNA (Fig. 3) 
and having identified MDA5 as the primary sensor involved in the detection of IBV (Fig. 4), we 
investigated whether there would be a temporal difference between IBV-induced accumulation of 
dsRNA and the upregulation of Ifnp transcription in CEK cells. Indeed, dsRNA could clearly be 
detected, even at low MOI of 0.01, by 12 hpi (Fig. 5A). In contrast, Ifnp levels at this time point 
remained low (Fig. 5B) even in cell cultures infected at higher MOIs of 1 or 10 and despite the 
increased abundance of dsRNA. To further investigate the time lag between early accumulation of 
dsRNA and late Ifnp expression, we performed a time course analysis. Foci of dsRNA could be 
detected as early as 3 hours post-infection only in IBV-infected cells (Fig. 5C, inset 3hpi) indicating 
that dsRNA starts accumulating very early in IBV-infected cells but apparently only leads to late 
(>18 hpi) Ifnp transcription. 

Primary CEK cells consist of a heterogeneous mix of cell types. Even at high MOI, IBV M41 infects 
only ~70% of the cells, indicating that not all cells are permissive to IBV M41 infection. In order to 
assess whether the time lag between accumulation of dsRNA and Ifnp expression could be due to 
the induction of Ifnp in bystander rather than IBV-infected cells, we used RNA fluorescent in situ 
hybridisation to visualise Ifnp mRNA in IBV-infected CEK cell cultures (Fig. 5D). At 12 hpi and low 
MOI (0.1), with most cells showing clear foci of dsRNA, none of the IBV-infected cells displayed an 
accumulation of Ifnp mRNA. At 12 hpi and a higher MOI, a few cells stained positive for Ifnp mRNA 
and only later, at 24 hpi, did most IBV-infected cells also stain positive for Ifnp mRNA, the kinetics 
of which closely following that observed in Figure 5A. In all cases, detection of Ifnp mRNA was 
restricted to cells that contained dsRNA. Altogether our data shows that IBV-infected, but not 
adjacent uninfected cells, upregulate Ifnp transcription in response to IBV infection. The significant 
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time lag between accumulation of dsRNA and Ifn/3 transcription further suggests the presence of a 
mechanism adopted by IBV to circumvent the onset of an IFN response. 

Accessory proteins 3a and 3b regulate IFN transcription and protein production 

To investigate whether the accessory proteins of IBV might play a role in the observed delay in Ifn/3 
transcription, we infected CEK cells with IBV scAUG3ab and scAUG5ab null viruses and the parental 
Beau-R virus (scAUG viruses possess a scrambled AUG start codon resulting in transcription but not 
translation of either ORFs 3a and 3b or 5a and 5b (36, 37)). Infection with the scAUG3ab, but not a 
scAUG5ab null virus resulted by 24 hpi in increased upregulation of Ifn/3 expression (Fig. 6A). 
Indicating that either one, or a combination of, accessory proteins 3a and 3b play a role in down 
regulating Ifn/3 transcription. The difference in Ifn/3 transcription between the scAUG3ab and the 
parental (Beau-R) virus could not be ascribed to differences in kinetics of virus replication, as all 
viruses displayed similar growth kinetics until 24 hpi (Fig. 6B). To determine whether 3a, 3b or 
both accessory proteins are involved in the observed down regulation of the IFN response, we 
quantified Ifn/3 transcription and IFN protein production in CEK cells infected with scAUG3a, 
scAUG3b and scAUG3ab mutant viruses, and compared the values observed in cells infected with 
Beau-R (Fig. 6C and 6D). Infection with all mutant viruses led to an increased transcription of Ifn/3 
when compared to the Beau-R (Fig. 6C), indicating that the presence of either one of the two 
accessory proteins is sufficient to limit Ifn/3 transcription. The kinetics of Ifn/3 transcription in 
response to AUG3a/b differs between Fig. 6A and Fig. 6C. In Fig 6C there is a significant difference 
in Ifn/3 transcription between AUG3a/b and Beau-R, which is absent in Fig 6A. This difference can 
probably be attributed to variation in the kinetics of Ifn/3 transcription between primary CEK cells 
isolated from embryos originating from different flocks. Nonetheless, this difference does not affect 
the conclusion that knockout of 3a and 3b leads to an increase in transcription of Ifn/3. 
No significant differences in IFN protein production were observed between cells infected with the 
Beau-R and the scAUG3ab double null virus, except at 36 hpi. Flowever, infection with scAUG3b 
virus led to an increase in IFN protein levels, whereas infection with the scAUG3a virus led to a 
decrease in IFN when compared to both Beau-R and scAUG3ab double null virus (Fig. 6D). Taken 
together, these results indicate that accessory proteins 3a and 3b both play a role in the inhibition 
of Ifn/3 transcription but have distinct and opposing effects on protein production. Accessory protein 
3b seems to be involved in limiting IFN protein activity whereas 3a is involved in promoting it. 
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Signalling of non-self dsRNA remains intact in IBV-infected cells 

Since IBV showed the intrinsic ability to delay Ifn/3 transcription in several cell types (Fig. 1 and 
Fig. 2) even in the presence of high levels of intracellular dsRNA (Fig. 5), we investigated the 
ability of IBV to interfere with sensing of non-self dsRNA by TLR3 or MDA5. We infected CEK cells 
with IBV M41 and subsequently used extracellular poly I:C to trigger TLR3 signalling (Fig. 7A). 
Stimulation with pI:C alone led to a significant increase in Ifn/3 transcription whereas stimulation 
with pI:C following an infection with IBV led to an enhanced increase in Ifn/3 transcription in an 
MOI-dependent manner. These results indicated that IBV infection does not interfere with TLR3- 
mediated Ifn/3 transcription, on the contrary IBV infection appears to result in a synergistic 
activation of the TLR3 pathway triggered by pI:C. Next, we investigated whether IBV infection 
could interfere with MDA5-mediated transcription of Ifn/3. Although transfection of pI:C into the 
intracellular compartment is a commonly used ligand of MDA5, this method induced very little 
transcription of Ifn/3 in primary CEK cells, because of low transfection efficiency (data not shown). 
As an alternative route to stimulate MDA5 in primary chicken cells, we investigated the use of 
either RVFV Cl 13 or IPNV, that induce Ifn/3 transcription in CEK cells (Fig 3C). RVFV Cl 13 is a (-) 
ssRNA virus with a truncated IFN antagonist (53), for which RIG-I, but not by MDA5 or TLR3, was 
previously shown to be the most likely PRR in mammalian cells (54, 55). Since chickens, as 
opposed to most mammals, do not have a RIG-I homologue, the most likely PRR for RVFV in CEK 
cells would be MDA5. IPNV is a birnavirus with a dsRNA genome that naturally infects salmonids 
but has been shown to enter but not replicate in cells of warm-blooded animals (56). To date, the 
PRR responsible for sensing IPNV dsRNA has not been described. Knockdown experiments in DF-1 
I7ti/3-Iuc reporter cells, using siRNAs against chicken MDA5 or TLR3, revealed that MDA5, but not 
TLR3, is the prime PRR for IPNV (Fig. 7B). These findings were confirmed using MEFs (mouse 
embryo fibroblasts) from knockout mice deficient in expression of either MDA5, RIG-I or the 
downstream adaptor protein MAVS. Here, knockout of either MDA5 or MAVS abrogated sensing of 
IPNV as shown by a strong reduction of Ifn/3 transcription, whereas knockout of RIG-I did not (Fig. 
7C). Both IPNV and RVFV Cl 13 were subsequently used to investigate whether IBV infection could 
interfere with MDA5-mediated transcription of Ifn/3 in CEK cells. 

Using quantification of Ifn/3 transcription by RT-qPCR as read out, we could show that IBV infection 
does not interfere with MDA5-mediated signalling of IPNV (Fig. 7D) or RVFV Cl 13 (Fig. 7E), in fact 
it had a synergistic effect on Ifn/3 transcription as previously observed for TLR3-mediated signalling 
(Fig. 7A). Similar results were obtained when stimulating IBV-infected DF-1 cells with IPNV or 


392 

393 

394 

395 

396 

397 

398 

399 

400 

401 

402 

403 

404 

405 

406 

407 

408 

409 

410 

411 

412 

413 

414 

415 

416 

417 

418 

419 

420 

421 

422 

423 


t[pI:C] (Fig. 7F), indicating that the observed synergistic effect is not specific to CEK cells. Taken 
together, IBV infection very efficiently prevents sensing of IBV (ds)RNA, but our results indicate 
that it does not interfere with sensing and downstream signalling of other non-self (ds)RNA ligands. 


DISCUSSION 

In this study we performed a comprehensive analysis of the kinetics of IBV infection in avian cells 
and studied the mechanisms by which IBV interferes with the onset of the type I IFN response. We 
show that infection with the Gammacoronavirus IBV leads to a considerable activation of the type I 
IFN response, albeit delayed with respect to the peak of viral replication and accumulation of viral 
dsRNA. Using an siRNA knockdown approach we show that MDA5 is the main receptor involved in 
the induction of Ifn/3 expression during IBV infection. We present evidence that IBV accessory 
proteins 3a and 3b play a role in the modulation of the delayed IFN response, by regulating 
interferon production both at the transcriptional as well as translational level. In addition, we show 
that although IBV alone effectively prevents Ifn/3 induction in IBV-infected cells, it does not block 
Ifn/3 induction upon stimulation of IBV-infected cells with other RIG-I, MDA5 or TLR3 ligands. To 
our knowledge, this study provides the most comprehensive analysis of the interplay between a 
Gammacoronavirus and the avian type I IFN response. 

Much of our current knowledge about the interaction of coronaviruses with the innate immune 
response (reviewed in (57)) comes from studies in mice and mouse cells using mouse hepatitis 
virus (MFIV). MFIV activated IFN production only in specific cell types and an efficient IFN response 
was only mounted by plasmacytoid dendritic cells (58), bone marrow derived macrophages (10, 
59) and oligodendrocytes (10). In a recent study on SARS-CoV and MERS-CoV in an epithelial lung 
cell line, ISGs started to be upregulated at 12 hpi (60), when virus titres were already reaching 
their maximum. The kinetics of IFN response observed in our study are in line with aforementioned 
studies, however it must be noted that in most cell types, infection with Alpha or Betacoronaviruses 
induced very little, if any, Ifn/3 transcription (8-13, 49). This suggests that all coronaviruses are 
able to modulate the activation of the type I IFN response. 

We found that IBV infection is detected by various chicken cell types, but until now it was unknown 
which PRR was involved. MFIV has been shown to be detected by MDA5 and not RIG-I or TLR3 in 
brain macrophages (50), by both MDA5 and RIG-I in an oligodendrocyte derived cell line (10) and 
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by TLR7 in plasmacytoid dendritic cells (58). Analysis of the chicken genome suggests that chicken 
lack a RIG-I homologue (52), and basal expression of Tlr7 was found to be very low in CEK cells 
(data not shown). We therefore silenced the remaining candidate RNA sensors MDA5 and TLR3, 
and were able to show that MDA5, but not TLR3, is involved in the sensing of IBV. Silencing of 
Mda5 did not lead to an increase in replication of IBV, suggesting that IBV might have developed 
strategies to counteract the activated IFN response. 

We recently reported membrane rearrangements in chicken cells infected with IBV (48), similar to 
those found in cells infected with Betacoronaviruses. In theory, the formation of intracellular 
membrane rearrangements might partly explain the discrepancy observed in the kinetics of dsRNA 
accumulation and Ifn/3 upregulation. Indeed, for SARS-CoV it has been shown that virus-induced 
double membrane vesicles (DMVs) contain dsRNA (14), suggesting that coronaviruses might 
exploit membrane structures to shield dsRNA from recognition by host PRRs (61). Flowever, the 
kinetics of Ifn/3 transcription were not investigated in these studies. The presence of coronavirus- 
induced DMVs has been demonstrated as early as 2 hpi in SARS-CoV-infected cells (14). Although 
we did not demonstrate the presence of DMVs in IBV-infected chicken cells at time points earlier 
than 7 hpi (48), it is likely that DMVs could also be present at earlier time points. As such, the 
timing of DMV formation in coronavirus-infected cells could suggest that membrane 
rearrangements play a role in the delayed activation of the IFN response by shielding dsRNA from 
cellular PRRs. 

In addition to membrane rearrangements, coronavirus-encoded proteins, including numerous 
accessory genes, have been shown to interfere with the type I IFN response pathway (reviewed in 
(17, 18). To investigate the possible role of IBV accessory proteins in the regulation of the IFN 
response we made use of our previously constructed mutant IBV Beau-R viruses that do not 
express either one or more of the four accessory proteins 3a, 3b, 5a and 5b. Previously we have 
demonstrated the accessory genes of IBV are not essential for replication (36, 37). In the present 
study we show that infection of CEK cells with 3a or 3b null viruses as well as a 3a/3b double null 
virus led to increased Ifn/3 transcription compared to Beau-R. Because the kinetics of Ifn/3 
transcription of 3a, 3b and 3a/3b null viruses are comparable to the parental virus, we conclude 
that, 3a and 3b are probably not responsible for the delay in Ifn/3 transcription, suggesting that IBV 
utilises additional strategies to delay transcription of Ifn/3. Apart from their effect on Ifn/3 
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transcription, 3a and 3b seem to have opposing effects on IFN protein production by IBV infected 
cells. Infection with the 3b null virus resulted in increased IFN production whereas infection with 
the 3a null virus resulted in reduced IFN levels compared to the Beau-R virus. Together with the 
observation that IFN production induced by the 3a/3b double null virus is comparable to that 
induced by Beau-R virus, our data suggests that accessory proteins 3a and 3b antagonise each 
other to tightly regulate IFN production (Fig. 6B). 

Using the eukaryotic linear motif server (62), we identified a Protein phosphatase 1 (PPl)-binding 
17 KISF 20 domain in the IBV 3b protein sequence. The canonical PPl-binding motif is 
[R/K][V/I/L]X[F/W], in which x can be any amino acid except proline (63). Interestingly, 
Alphacoronavirus TGEV accessory protein 7 (TGEV-7) has been shown to bind PP1 via a binding 
motif similar to that found in IBV 3b (64). Similar to IBV scAUG3b, infection with TGEV-A7 led to 
increased mRNA and protein levels of IFN(3 (65). The fact that both TGEV-7 and IBV 3b contain a 
PP1 binding domain indicates that interaction with PP1 could be a common strategy of 
coronaviruses to inhibit the host innate immune response. The mechanism by which interaction of 
coronavirus accessory proteins with PP1 counteracts the innate immune response still needs to be 
determined. One clue might come from the PPl-binding domain of Measles virus V, which was 
recently shown to be essential for inhibition of MDA5 signalling (66, 67). Measles V protein binds 
PP1 and inhibits dephosphorylation of MDA5, which is required for activation and subsequent 
signalling by MDA5. Motif analysis for IBV 3a protein did not reveal the presence of relevant 
motives that might explain the observed activity of 3a on IFN regulation. We conclude that both 
accessory proteins 3a and 3b limit Ifn/S transcription but have distinct and opposing effects on 
protein production. Whilst 3a seems to promote IFN production, 3b seems to be involved in limiting 
IFN protein production, possibly through a similar mechanism as described for protein 7 of TGEV. 
The fact that IBV 3a and 3b have opposing roles in regulating IFN production indicates that CoV's 
tightly regulate IFN production to balance their own survival with that of the host. This hypothesis 
is supported by the observation that field isolates lacking 3a and 3b display reduced virulence in 
vitro as well as in vivo (68). Elucidation of the exact mechanisms of action of 3a and 3b will be the 
subject of further investigation. 
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To investigate whether IBV interferes with a general sensing of (ds)RNA ligands or downstream 
signalling that leads to Ifn/3 transcription, we stimulated IBV-infected cells with TLR3, RIG-I, and 
MDA5 ligands. Surprisingly, we found that infection with IBV did not reduce Ifn/3 transcription but 
rather increased Ifn/3 levels upon stimulation with these PRR ligands. Similar to IBV, MHV has been 
shown unable to inhibit expression of Ifn/3 induced by either t[pI:C] or Sendai virus (69, 70), but in 
these studies no synergistic effect was observed. Currently, we can only speculate about the cause 
of this synergistic effect. It appears that IBV infection 'arms' the Ifn/3 induction pathway, without 
actually triggering it, possibly by enhancing the activity of one or more components of the pathway 
leading to Ifn(3 upregulation. One possibility is that IBV-proteins interact with host-proteins that 
regulate this pathway through ubiquitination and phosphorylation (reviewed in (71)). The fact that 
stimulation with either TLR3 or MDA5 ligands resulted in exacerbated transcription of Ifn/3 indicates 
that IBV influences a component which is downstream of both MDA5 and TLR3. 

Taken together, our study provides the first comprehensive analysis of host-virus interactions of a 
Gammacoronavirus with the avian innate immune response. We show that the Gammacoronavirus 
IBV, induces activation of the type I IFN response in primary chicken renal cells, tracheal epithelial 
cells and in a chicken cell line. We show that activation of the IFN response is dependent on MDA5 
but is delayed with respect to the peak of virus replication. We demonstrate that Ifn/3 transcription 
is restricted to IBV-infected, dsRNA-containing cells and provide evidence that accessory proteins 
3a and 3b of IBV are involved in regulating transcription as well as protein production of type I IFN. 
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FIGURE LEGENDS 


FIG 1. IBV infection delays Ifn/3 upregulation. 

Chicken embryo kidney (CEK) cells were infected with IBV M41 at an MOI of 0.1. (A) Replication of 
IBV was quantified by titration in cell culture supernatants of infected cells; in a parallel 
experiment, intracellular IBV RNA was quantified using RT-qPCR. (B) 17/7/3 mRNA levels were 
determined using RT-qPCR and IFN protein levels using a chicken IFN-specific mx-luc cell-based 
bioassay, respectively. (C) Expression of genes involved in the antiviral response. All gene 
expressions were calculated as fold changes relative to uninfected control cells and normalised 
against an external reference gene (luciferase). For IBV total RNA, fold changes were calculated 
relative to C t 30. Depicted are the results of a representative experiment out of three independent 
experiments. 

FIG 2. Delayed induction of Ifn/3 transcription is independent of cell type or IBV strain. 

(A) Epithelial cells from adult chicken trachea were infected with IBV M41 at MOI 0.1. (B) 
Fibroblast DF-1 cells were infected with IBV Beau-R at MOI 0.1. (C) CEK cells were infected with 
IBV strains M41, QX and It02, at an MOI 0.1. Intracellular IBV total RNA (open diamonds) and Ifn/3 
mRNA (bars) are depicted as fold changes as assessed by RT-qPCR. Gene expression of Ifn/3 was 
calculated as fold changes relative to uninfected control cells and normalised against an external 
reference gene (luciferase).. For IBV total RNA, fold changes were calculated relative to C t 30. Error 
bars indicate standard deviation. 

FIG 3. Chicken cells have the intrinsic ability to respond rapidly to dsRNA. 

(A) CEK cells were seeded in 24 well plates and 48 hours later stimulated with extracellular poly 
I:C for the indicated times. (B) DF-1 cells were infected with IPNV, a non-replicating dsRNA virus, 
or stimulated with extracellular pI:C or transfected pI:C (t[pl:C]). Four hours later, Ifn/3 fold 
changes were determined by RT-qPCR. Bars represent the mean (plus standard deviation) of 
triplicate wells from a representative experiment. Asterisks indicate significant differences (P<0.01) 
with respect to the non-stimulated control as assessed by one-way ANOVA followed by a Bonferroni 
post-hoc test. (C) CEK cells were infected with IBV M41 (MOI 1), IBV Beau-R, (MOI 1), Sindbis-GFP 
(MOI 1), IPNV (MOI 50) and RVFV Cl 13 (MOI 5). Depicted are Ifn/3 fold changes at 12 hpi relative 
to uninfected control cells as assessed by RT-qPCR. 
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FIG 4. MDA5, and not TLR3, is the prime sensor of IBV. 

(A) CEK cells were infected with IBV M41 for 24 hours, in the presence or absence of RNase A. Ifn/3 
expression was analysed by RT-qPCR. Stimulation with pI:C in the presence or absence of RNase A 
was included as a positive control. (B,D) DF-1 cells and (C) DF-1 Ifnp-luc reporter cells were 
transfected with siRNAs against Tlr3, Mda5 or a control siRNA and 48 hours later infected with IBV 
M41 (MOI 0.1). (B) Ifn/3 mRNA, (C) Ifnpi-luciferase activity, and (D) IBV titres and intracellular RNA 
were analysed 18 hpi. Bars represent the mean (plus standard deviation) of triplicate wells from a 
representative experiment. Asterisks indicate significant differences (P<0.01) with respect to the 
non-RNaseA-treated control (A) or to the siRNA control (B-C), as assessed by one-way ANOVA 
followed by a Bonferroni post-hoc test. 

FIG 5. Early accumulation of dsRNA in IBV-infected cells does not result in early 
induction of Ifn/3 

CEK cells were infected with IBV M41 or IBV Beau-R at the indicated MOIs. At time point 12 hpi (A) 
dsRNA was visualised in M41-infected cells using an antibody against dsRNA. (B) Expression of Ifn/3 
mRNA was analysed by RT-qPCR. (C) CEK cells were infected with IBV M41 and accumulation of 
dsRNA was visualised at the indicated time post infection. (D) RNA fluorescent in situ hybridisation 
of Ifn/3 mRNA in IBV M41-infected CEK cells. Open arrowheads indicate cells that contain dsRNA 
and no Ifn(3 mRNA. Solid white arrowheads indicate cells that contain both dsRNA and Ifn/3 mRNA. 

FIG 6. Accessory proteins 3a and 3b are involved in regulation of IFN transcription and 
protein production. 

(A) CEK cells were infected with IBV Beau-R 3a/3b (scAUG3ab) or 5a/5b (scAUG5ab) null viruses 
(MOI 0.1). Ifn/3 levels were determined using RT-qPCR. (B-D) CEK cells were infected with 
scAUG3a, scAUG3b or scAUG3ab null IBV viruses (MOI 0.1). In the same cultures (B) Ifn/3 mRNA, 
(C) virus titres and (D) type I IFN protein were quantified. Bars represent the mean (plus standard 
deviation) of triplicate wells from a representative experiment. Significant differences (P<0.01) 
relative to the Beau-R virus at the same timepoint (*) or between the indicated bars (#) as 
assessed by a two-way ANOVA followed by a Bonferroni post-hoc test. 
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FIG 7. Signalling of non-self RNA remains intact in IBV-infected cells 

(A) CEK cells were infected with IBV M41 for 3 hours and stimulated with extracellular poly I:C (50 
pg/ml) for an additional 3 hours after which Ifn/3 transcription was analysed by RT-qPCR. (B) DF-1 
IfnfHuc reporter cells were transfected with siRNAs against Tlr3, Mda5 or a control siRNA and 48 
hours later infected with IPNV (MOI 50); at 6 hpi luciferase activity was quantified. (C) Knockout 
(KO) and wild-type (wt) MEFs were infected with IPNV (MOI 50) for 8 hours. (D) CEK cells were 
infected with IBV M41 (MOI 10) for 6h and super-infected with IPNV or UV-inactivated IPNV (MOI 
50) for an additional 6 h. (E) CEK cells were co-infected with IBV M41 (MOI 5) and RVFV clone 13 
(MOI 5) and sampled at 6 hpi. (F) DF-1 cells were infected with IBV Beau-R (MOI 1) for 3 h and 
super-infected with IPNV (MOI 50) or transfected with pI:C (t[pI:C], 500 ng/well) for an additional 
4 h. (C-F) Ifnfi levels were quantified by RT-qPCR. Bars represent the mean (plus standard 
deviation) of triplicate wells. Significant differences (P<0.01) are indicated by (*) as assessed by 
one-way ANOVA followed by a Bonferroni post-hoc test. 


gene 


kind 


sense sequence (5'-3') 


Accession nr. 


reference 


TLR3 

siRNA 

S 

AS 

UCGAAUACUUGGCUUUAAA 

UUUAAAGCCAAGUAUUCGA 

NM_001011691 


Ctrl 

siRNA 

S 

AGGUAGUGUAAUCGCCUUG 

- 




AS 

CAAGGCGAUUACACUACCU 



MDA5 

siRNA 

S 

ACACUGGUAUCAAGUUAUU 

GU570144 




AS 

AAUAACUUGAUACCAGUGU 



IFN[3 

RQ primer 

FW 

GCTCTCACCACCACCTTCTC 

ENSGALT00000039477 




RV 

GCTTGCTTCTTGTCCTTGCT 



IFNa 

RQ primer 

FW 

ATCCTG CTG CTC ACG CT CCTT CT 

XM_004937096 

40 



RV 

GGTGTTGCTGGTGTCCAGGATG 



IRF3 

RQ primer 

FW 

CAGTGCTTCTCCAGCACAAA 

NM_205372 




RV 

TGCATGTGGTATTGCTCGAT 



IRF1 

RQ primer 

FW 

CAGGAAGTGGAGGTGGAGAA 

ENSGALG00000006785 




RV 

TGGTAGATGTCGTTGGTGCT 



TLR7 

RQ primer 

FW 

TTCTGGCCACAGATGTGACC 

NM_001011688 

40 



RV 

CCTT CAACTT GGCAGTGCAG 



TLR3 

RQ primer 

FW 

TCAGTACATTTGTAACACCCCGCC 

NM_001011691 

40 



RV 

GGCGT CATAAT CAAACACT CC 



MDA5 

RQ primer 

FW 

TGGAGCTGGGCATCTTTCAG 

GU570144 




RV 

GTT CCCACGACT CT CAATAACAGT 



Mx 

RQ primer 

FW 

TTGTCTGGTGTTGCT CTT CCT 

ENSGALT00000025999 




RV 

GCTGTATTTCTGTGTTGCGGTA 



OAS 

RQ primer 

FW 

CACGGCCTCTTCTACGACA 

NM_205041 

41 



RV 

TGGGCCATACGGTGTAGACT 



IL8 

RQ primer 

FW 

TTGGAAGCCACTTCAGTCAGAC 

NM_205498 

41 



RV 

GGAGCAGGAGGAATTACCAGTT 



PKR 

RQ primer 

FW 

CCTCTGCTGGCCTTACTGTCA 

NM_204487 

42 



RV 

AAGAGAGGCAGAAGGAATAATTTGCC 


ADAR 

RQ primer 

FW 

TGTTTGTGATGGCTGTTGAG 

AF403114 




RV 

AGATGTGAAGTCCGTGTTG 



ISG12 

RQ primer 

FW 

TAAGGGATGGATGGCGAAG 

NM_001002856 




RV 

GCAGTATCTTTATTGTTCTCAC 



MHC-I 

RQ primer 

FW 

CTT CATT GCCTTCGACAAAG 

NM_001031338 

41 



RV 

GCCACTCCACGCAGGT 



IFNAR2 

RQ primer 

FW 

GCTTGTGTTCGTCAGCATT 

ENSGALT00000036778 

41 



RV 

TT CGCAAT CTT CCAGTT GT 



IBV-N 

RQ primer 

FW 

GAAGAAAACCAGTCCCAGA 

AY851295 




RV 

TTACCAGCAACCCACAC 



Luciferase RQ primer 

FW 

TGTTGGGCGCGTTATTTATC 

X65316 




RV 

AGGCTGCGAAATGTTCATACT 
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